Electronic transport in the Coulomb phase of the pyrochlore spin ice 
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We investigate the transport properties of itinerant electrons interacting with a background of 
localized spins in a correlated paramagnetic phase of the pyrochlore lattice. We find a residual resis- 
tivity at zero temperature due to the scattering of electrons by the static dipolar spin-spin correlation 
that characterizes the metallic Coulomb phase. As temperature increases, thermally excited topo- 
logical defects, also known as magnetic monopoles, reduce the spin correlation, hence suppressing 
electron scattering. Combined with the usual scattering processes in metals at higher temperatures, 
this mechanism yields a non-monotonic resistivity, displaying a minimum at temperature scales 
associated with the magnetic monopole excitation energy. Our calculations agree quantitatively 
with resistivity measurements in Nd2lr207, shedding light on the origin of the resistivity minimum 
observed in metallic spin-ice compounds. 



The interaction between itinerant electrons and lo- 
calized moments continues to attract considerable in- 
terest as model systems to understand non-Fermi liq- 
uid behaviors and unusual transport phenomena. Since 
the pioneering work of Kondo on the resistivity min- 
imum in metals [ij, the Kondo-lattice model has be- 
come the paradigm for a new class of problems [H . This 
simple Hamiltonian describing the competition between 
Kondo screening of individual moments and direct ex- 
change among localized spins plays an important role 
in the physics of heavy- fermion materials 0]. Novel 
magneto-transport phenomena such as colossal magne- 
toresistance [J and anomalous Hall effect [HI emerge in 
a closely related double-exchange model in which the 
Kondo-screening is replaced by the effect of a ferromag- 
netic Hund's rule coupling 

Adding geometrical frustration to Kondo-lattice mod- 
els opens a new avenue for exploring unconventional 
phase transitions and many-body states. For example, 
a spin-liquid phase with broken time-reversal symmetry 
and spontaneous Hall effect was recently observed in the 
pyrochlore iridate Pr2lr207 0- The conduction 5d elec- 
trons in this compound coexist with the localized 4/ mo- 
ments at the Pr sites that form a frustrated pyrochlore 
lattice (Fig. [T]). A strong easy-axis anisotropy forces the 
moments to point along the local (f 11) directions. When 
spin interactions are restricted to nearest-neighbor fer- 
romagnetic exchange Jp, an extensively large number of 
Ising configurations satisfying the "ice rules" are degener- 
ate ground states [§]. As temperature tends to zero, this 
degeneracy is eventually lifted by residual perturbations, 
such as spin-electron coupling [10|. However, spatial spin 
correlations exhibit an unusual power-law decay over a 
wide temperature range below the Curie- Weiss constant 
©cw ~ Jp/ks- Elementary excitations in this so-called 
Coulomb phase [llj are defect tetrahedra that violate the 



ice rules and carry a finite magnetic charge |12. 13|. 

Recently, an unexpected resistivity minimum, similar 
to the one seen in Kondo systems, was observed in the 
paramagnetic phase of the metallic pyrochlore iridates 
Pr2lr207 and Nd2lr207 [lj,ll5|. Although quantum fluc- 



tuations of these Ising-like moments were shown to be 
non-negligible jTij . the large easy-axis anisotropy sup- 
presses the Kondo screening. On the other hand, a non- 
Kondo mechanism for the resistivity minimum was sug- 
gested by a recent numerical study of a strong-coupling 
Kondo-lattice model assuming classical Ising spins on the 
pyrochlore lattice 17|. In particular, their cellular dy- 
namical mean-field calculation suggests a clear correla- 
tion between the position of resistivity minimum and the 
onset temperature of ice rules. However, up to date, there 
is not a transparent theory relating the resistivity of the 
Coulomb phase to its unique magnetic correlations and 
excitations. 

In this paper, we present an analytical theory for elec- 
tron transport in the Coulomb phase of spin ice. In the 
continuum approximation, the problem is recast into that 
of a correlated random "magnetic" field locally coupled 
to a Fermi sea of electrons. The resulting electronic scat- 
tering is the dominant transport mechanism at low tem- 
peratures. We show that thermally excited monopoles 
reduce the spin-spin correlation and suppress the elec- 
tron scattering rate. Combined with the other scattering 
processes taking place in the metallic compounds - due, 
for instance, to electron-phonon interaction - we obtain 
a resistivity minimum at temperature scales associated 
with the excitation energy of the magnetic monopoles. 
We compare our calculations with the resistivity data on 
Nd2lr207 and find a good quantitative agreement. Our 
results show that strong geometrical frustration com- 
bined with topological excitations give rise to a resistivity 
minimum of non-Kondo character in metallic spin ice. 
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FIG. 1: (Color online) A fragment of the pyrochlore lattice. 
In the Coulomb phase, local spin arrangement on each indi- 
vidual tetrahedra obeys the so-called ice rules with two spins 
pointing in and two spins pointing out. 



We start with a ferromagnetic Kondo-lattice model on 
the pyrochlore lattice: 



,Cj-„ +h.c. 



-Jh E ^» ■ c|,aCi,/3 - S, • Sj. (1) 



Here the first term describes electron hopping between 
nearest-neighbor sites, t is the hopping integral, and c| ^ 
creates an electron with spin a , J, on site i. The itiner- 
ant electrons interact with the localized spins S,; through 
an on-site Hund's coupling Jh; here aajs is a vector of 
Pauli matrices. As discussed above, the magnetic mo- 
ments Si = TiS Bi are forced to point along the local 
Gi = (111) axes by a strong easy-axis anisotropy; n = ±1 
is an Ising variable. The last term represents a ferromag- 
netic {Jp > 0) exchange interaction between neighboring 
moments. 

We first consider the ground state of the localized spin 
system. The exchange energy of a single tetrahedron is 
minimized by six different Ising states with two spins 
pointing in and two pointing out of the tetrahedron. 
This 2-in-2-out rule is analogous to the Bernal-Fowler 
rule for water ice [§]. As explained by Pauling, the very 
large number of configurations satisfying the ice rule in 
a macroscopic sample gives rise to a measurable resid- 
ual entropy at low temperatures Because of this 
huge degeneracy, spins remain disordered even at tem- 
peratures well below the exchange energy scale Jp. 

Remarkably, spatial correlation of spins in this disor- 
dered yet highly constrained phase exhibits a power-law 
decay. This is because the local 2-in-2-out ice rule trans- 
lates into a divergence-free condition V • B = for the 
coarse-grained magnetization field B(r) = J2iev ^i/^i 
where Vr denotes a coarse-graining block of volume 



n pp'^ that contains several tetrahedra. In momen- 
tum space this constraint becomes k • B(k) = 0, indicat- 
ing that only transverse fluctuations are allowed. Con- 
sequently, the static spin correlator has a dipolar form 
(B^(r)S,(0)) - (Sf,, - 3f^f,)/r3 for large r. 

The low-energy excitations of this paramagnetic phase 
are topological defects that carry magnetic charges, 
namely, tetrahedra with 3-in-l-out or l-in-3-out spin con- 
figurations. Since these emergent magnetic monopoles 
are sources and sinks for the magnetization field, hop- 
pings of thermally excited monopoles modify the spin 
correlator at finite temperatures. In the low-T hydro- 
dynamic regime, the relaxation of the magnetization by 
monopole motion is governed by a stochastic differential 
equation [11, ^ 



dB 

'dt 
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Here n„i = V • B is the monopole density, /i and D are 
the monopole mobility and diffusion constant, respec- 
tively, C{t) is a Gaussian noise source, $ = (16/y/3) aT 
indicates the entropic origin of monopole drift |12| . and 
a is the lattice constant. The first term in the right- 
hand side of Eq. ^ implies a magnetic relaxation time 
Tm ^ ^/fJ'n-m^- Since the mobility ^ ~ 1/T at low tem- 
peratures [igj . the combination ^$ is a relatively weak 
function of T, and the relaxation time is governed by the 
monopole density: ~ e~^l'^ and diverges exponen- 
tially Tm(T) = T,„ e^^"^ HSl; here A is the activation 
energy for creating monopoles. Eg. (|^ gives rise to a 
time-dependent correlation function |lll. |21|: 
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where K is the stiffness constant of the fiux field, 
= Id oc jrim is the Debye screening length of 
monopoles. The two terms in Eq. ([3]) denote the trans- 
verse and longitudinal components, respectively. 

Since the transport properties are dominated by con- 
duction electrons near the Fermi surface, their wavcfunc- 
tion has a characteristic wave length •p'^ which is much 
larger than the underlying lattice constant. Considering 
the low electron density of the iridate compounds, we can 
choose a coarse-graining block V whose volume ~ P~p' ■ 
We thus consider a model Hamiltonian describing elec- 
trons interacting with the averaged magnetization field: 

•H = ^ £p 4_„Cp - |: E ^('^^ ■ c\ o,C^+^,,fi■, (4) 



p,a 



p,k,a/3 



where is the electron energy measured with respect 
to the chemical potential ^l, g Jh is the coupling 
constant, and the correlation for the magnetization field 
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is given by Eq. ([3]). For simplicity, here we consider a 
parabolic dispersion Sp = /2m — /i, where m is the 
effective electron mass. Note that the Hamiltonian is 
weakly time-dependent owing to the slowly varying mag- 
netic field. We first calculate the lifetime of electrons 
close to the Fermi surface. To this end, we consider 
the leading second-order terms of the electronic self- 
energy corresponding to the diagram shown in Fig. [2] 
Here the solid line denotes the free electron propagator 
G'o(p,w) = {oj — Ep ± iO)~^ with ± sign corresponding 
to advanced/retarded Green's function, respectively, and 
the wavy line represents the B-ficld correlator. Sepa- 
rating the correlator into the transverse and longitudi- 
nal components, we obtain a self-energy Sa.Q/3 — ^ap'^a 
whose imaginary part is given by 

ImS,(p,u;) = -4^ V / dkkVa{k) 

X tan^""^ {Ta[vk — s{u} ~ Ep ~ fc^/27n)]} , (5) 

Here a = T, L labels the transverse and longitudinal 
components of the correlator, v = p/m is the electron 
velocity, A = is an ultraviolet cutoff, Vt = 2 and 
2?L = /{k'^ + K^). The magnetic relaxation times arc 
tt — Tm and tl — 4- D|kp). Since the cor- 

relation (jSj is valid only for coarse-grained magnetiza- 
tions, the length scale should be over several lattice 
constants. A natural choice for the cutoff is the Fermi 
wavevector pp. We note that our results do not depend 
on the exact value of the cutoff, as we show below. 

We now consider the decay rate t^^^ = 2ImI] of elec- 
trons on the Fermi surface w = ep = 0. Fig. [3] shows 
the transverse and longitudinal components of ImE as a 
function of temperature. In particular, the integral ^ 
for the transverse part depends only on two dimension- 
less parameters A/pp and ep Tm,Oj where ep is the Fermi 
energy. The longitudinal part, on the other hand, de- 
pends also on the monopole screening length ijj. How- 
ever, since at low temperature £p) ^ ic and 2?^ <C 
for most of the integration region, the longitudinal con- 
tribution to the decay rate ImE is negligible compared 
with the transverse part. 

Interestingly, the transverse part ImEr decreases with 
increasing T. This is in stark contrast to the behavior 
of the scattering rate due to other common processes in 
metals, such as scattering by phonons, electrons, or spin 
fluctuations. The anomalous temperature dependence of 

Go(p,a;) \ 



FIG. 2: Lowest-order diagram of the electron self-energy 
"Ea^ {p,u). The solid and wavy lines denote the electron prop- 
agator Go and the spin correlator @, respectively. 
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FIG. 3: (Color online) (a) Transverse and (b) longitudinal 
components of the electron scattering rate Te~^(T) — 2ImE as 
a function of temperatue. Both are normalized with respect to 
ImET(r = 0) = g^l\^ /2-kvfK. We have used Tm,oDp\. = 0.1 
and a Debye screening length Id{T = A) = 1.65 x 10"^ in 
the calculations. 

ImS can be understood as follows. At T = 0, a constant 
relaxation time results from scattering of electrons by 
a dipolar-correlated static random field. At finite tem- 
peratures, thermally excited monopoles reduce the spa- 
tial spin correlation and suppress the electron scattering 
rate. In particular, expanding ImI]T(pF,w = 0) at low 
temperatures gives 

2TTVpK \ ttA^ epTmfl 2pp J 

(6) 

that decreases exponentially at small T. It should also 
be noted that at high temperatures T ^ A, the correla- 
tor ([3]) of the coarse-grained magnetization is no longer 
valid due to proliferation of monopoles. Instead, a power- 
law relaxation rate oc T" results from electron- 
phonon or electron-magnon scatterings. The competition 
of these two scattering mechanisms thus gives rise to a 
resistivity minimum. 

We now turn to the calculation of the resistivity in 
the Coulomb phase, using the so-called memory func- 
tion approach [i^. In this method the conductivity is 
expressed as a{z,T) — {iujp/4:7r)[z + M{z,T)]~^ , where 
ujp = Aire^ne/m is the electronic plasma frequency, Ue 
and m are the electron density and mass, respectively, 
and the memory function is approximated as M{z) = 
(m/riez) [(t>{z) — (/>(0)] with 

/>oo 

^{z) = -z dte^^'{[[%,j]u[H,j]t=o\). (7) 
Jo 

Here j = = Z]p ^v(p)Cp,QCp,a current operator 

(/i = X, y, z), w^(p) = dep/dpfi = Pfi/rn, the Hamiltonian 
H is given by Eq. ([4]), and (• • • ) denotes average over the 
electronic ground state and degenerate spin configura- 
tions with correlations specified by Eq. ([3]). 
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Since, as discussed above, the scattering is dominated 
by the transverse part of the spin correlation, we thus 
focus on its contribution to the conductivity in the fol- 
lowing discussion. A straightforward calculation gives 



M{z) 



SmUeK J (27r)3 (27r)3 ' ' 
/(gp+k) - /(gp) 



(8) 



{z + Ep- Ep+k + i/Tm) (£p+k - Ep - i/Tm) ' 

where /(ck) denotes the Fermi function of energy Sk- 
Fig. (U) shows the resulting dynamical conductivity as 
a function of frequency uj for various temperatures. By 
expanding M{z) at low temperatures in powers of lo and 
exp(— A/T), wc obtain the standard Drude Lorentzian 
form [H: cr(w) = a{ujp/ ATT){iuj + t^^^) / {uj"^ -^t^^"^), where 
a = (1 + 9M79w|„=o)~^ and = l/aM"(0) is the 
transport relaxation time. The dc resistivity is related 
to the transport lifetime as piT) = {in/auip) (T). 
Fig. shows the dc resistivity p{T) vs temperature 
for various dimensionless parameter rj = ep Tm,o- Similar 
to the electronic relaxation time, the resistivity increases 
to a maximum value as T tends to zero. At low temper- 
atures we again find that p{T) decreases exponentially 
with T: 



PiT) = po[l- 



8 pf e 



-A/T 



StT a EpTn 



(9) 



where po oc [g^ / K){p']pe p / UeUj'^) is the residual resistivity 
at zero temperature. 

To leading linear approximation, we use the 
Matthiessen's rule to combine various contributions to 
the resistivity. In particular, we consider a temperature 
dependence: ptot{T) ~ po - Ae-^/"^ + BT", where A is 
the coefficient of the exponential term in ([9]) and B > 0. 
For example, the exponent n = 2 for scattering due to 
electron-electron interaction in a Fermi liquid. We find 
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FIG. 4: (Color online) Dynamical conductivity a{u},T) as a 
function of frequency to at various temperatures. The con- 
ductivity is measured with respect to ao = uip/inA. 
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FIG. 5: (Color online) (a) Resistivity p{T) (in units of po) as 
a function of temperature T for varying dimensionless param- 
eter 77 = spTm.o- The Fermi energy is fixed at ef = lOOA 
in the calculation, (b) The rescaled resistivity r]{p — po) col- 
lapses on a universal curve at low temperatures. Also shown 
is the rescaled experimental data from the measurements on 
NdalraOr M- 



that the resistivity has a minimum at 

A/(n+f) 



W 



( ^A ( Bn\irvT\ 
\n+l \AA) ) 



(10) 



for A> Ac^ nS/A(A/(n-M)e)"+\ Here W{x) de- 
notes the Lambert M^-function. Since —1 < W[—x) < 
for a; > 0, the resistivity minimum thus occurs at tem- 
perature scales associated with the magnetic monopole 
excitation e nerg y, consistent with the numerical calcula- 
tion in Ref. [U . 

Experimentally, the resistivity upturn observed in 
Nd2lr207 under pressure agrees well with our calcula- 
tions. Eq. (|9]) shows that the low-T resistivity curves 
vip ^ Po) rescaled by the dimensionless parameter t] = 
£f Tm,,o collapse on a universal curve independent of the 
magnetic relaxation time. By plotting the rescaled ex- 
perimental resistance i?exp — Rq as a function of T/A, 
the data points indeed fall on the theoretical curve at 
low temperatures [Fig. EJb)]. From this analysis, we 
obtain a monopole activation energy A « 5.2 K, inde- 
pendent of applied pressures, and a residual resistance 
Rq « 0.87 $7 and 0.52 O for samples under pressure 6 
and 9 GPa, respectively. The ratio of the dimension- 
less parameters ry(P) at the two pressure values P is 
7? (6 GPa) /ry (9 GPa) — 3, i.e. the prefactor of the mag- 
netic relaxation rate l/r™ o increases by a factor of three 
as the pressure is increased from 6 to 9 GPa. This is 
not unexpected, since the spin-spin separation shrinks 
under pressure and, consequently, the exchange coupling 
increases. We also note that, while the resistivity mini- 
mum observed in the other compound Pr2lr207 is com- 
patible with our model prediction, the low temperature 
data seems to suggest a In T behavior [14{ , implying the 
effect of partial Kondo-screening. 
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In summary, we have shown that electronic trans- 
port in spin ice exhibits a resistivity minimum that docs 
not originate from Kondo screening. Instead, it comes 
from the scattering of electrons by a random magne- 
tization field with long-range correlation, which is con- 
trolled by thermally excited topological defects, or emer- 
gent magnetic monopoles. Using the memory function 
approach, we have shown that the dc resistivity decreases 
with increasing temperatures, in agreement with exper- 
iment and with numerical simulations of the Kondo- 
lattice model. Our theory opens a new route to explore 
unusual transport phenomena in metallic geometrically 
frustrated magnets. Questions that remain open, and 
deserve further study, are how the magnetic relaxation 
rate is affected by quantum spin dynamics and the feed- 
back effect of the itinerant electrons on the local moments 
correlations. 
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